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Interband optical transitions in graphene are subject to pseudospin selection rules. Impulsive 
excitation with linearly polarized light generates an anisotropic photocarrier occupation in momen¬ 
tum space that evolves at timescales shorter than lOOfs. Here, we investigate the evolution of 
non-equilibrium charges towards an isotropic distribution by means of fluence-dependent ultrafast 
spectroscopy and develop an analytical model able to quantify the isotropization process. In con¬ 
trast to conventional semiconductors, the isotropization is governed by optical phonon emission, 
rather than electron-electron scattering, which nevertheless contributes in shaping the anisotropic 
photocarrier occupation within the first few fs. 


Introduction — The unique optical and electronic prop¬ 
erties of graphene are appealing for advanced appli¬ 
cations in photonics and optoelectronics pj, [gj. A va¬ 
riety of prototype devices have already been demon¬ 
strated, such as transparent electrodes in displays [ij 
and photovoltaic modules^, optical modulatorsplas- 
monic devices (il-lioj. microcavities [Til IT^ and ultra¬ 
fast lasers [TiJ. Amongst these, a significant effort 
is being devoted to the development of broadband 
photodetectors [ 3 ]. In this context, the same properties 
that make this two-dimensional (2d) material so appeal¬ 
ing are also profoundly different from standard semicon¬ 
ductors. Thus, in order to fully exploit the technological 
potential of graphene, a full and comprehensive descrip¬ 
tion of the physical phenomena occurring when light ex¬ 
cites transitions within the Dirac cones is needed. 

The ultrafast carrier dynamics has been extensively in¬ 
vestigated in graphene 15- lfj by tracking the evolution 
of the non-equilibrium distribution created by impulsive 
optical excitation. Once a Fermi-Dirac distribution is 
established [3, Ho), it cools down by optical phonon emis¬ 
sion at the hundred fs timescale [16|, [2JJ and by scattering 
on acoustic phonons [22l - [24| at the ps timescale. How¬ 
ever, the complete description of the phenomena respon¬ 
sible for the momentum-space photocarrier redistribution 
within the first lOOfs requires significant efforts: a com¬ 
plex interplay between electron-electron (e-e) 18, 20| and 
electron-phonon (e -Ph)0,0 scattering dominates the 
thermalization dynamics. 

Several theoretical[25l |26l and experimental [27l l28j 
studies demonstrated the possibility to generate an 
anisotropic photocarrier distribution in the momentum 
space by means of linearly polarized light, and exploit¬ 
ing the pseudospin selection rules, which follow from 
the two sublattices present in graphene [3] . The pseu¬ 
dospin can be seen as an internal angular momentum j3pj . 
Thus, similarly to a gyroscope, a torque has to be ap¬ 
plied in order to change the pseudospin orientation[30|. 
This can be achieved by the interaction with an electro¬ 
magnetic wave whose electric field provides the necessary 



FIG. 1. (a) Anisotropic carrier distribution after excitation 
with linearly polarized light due to pseudospin-momentum 
coupling, (b) A linearly polarized pump pulse is followed, 
after time delay At, by a probe pulse with polarization ro¬ 
tated by 45°. After interaction with graphene, a polarizing 
beam splitter separates parallel (Em) and perpendicular (Ej_) 
components of the probe field for simultaneous detection. 


momentum [dll. The torque equals zero when the pseu¬ 
dospin is parallel to the electric field polarization making 
its flip forbidden. In contrast, this process is most ef¬ 
ficient when the light polarization and pseudospin are 
normal to each other [30j. and the torque is maximized. 
Thus, as shown in Fig Oja), a pseudospin flip is neces¬ 
sary to directly excite an electron from the valence to 
the conduction band. Hence, an excitation with linearly 
polarized light results in an anisotropic carrier occupa¬ 
tion within the Dirac cone. 

Here we measure the real-time isotropization of the 
photocarrier distribution in graphene for different con¬ 
centrations of the out-of-equilibrium carriers and explain 
the main mechanism responsible for the anisotropy re¬ 
laxation. We perform fluence-dependent measurements 
of the carrier anisotropy with few-fs temporal resolution 
and develop an analytical description of the transient op- 
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tical absorption. This allows us to resolve the pseudospin 
dynamics over its timescale. The combination of exper¬ 
iment and theory allows us to demonstrate that (i) the 
photocarrier distribution anisotropy is stronger for lower 
radiation intensities when the anisotropic character of the 
excitation is not yet suppressed by the isotropic Fermi- 
Dirac sea; (ii) the isotropization dynamics is driven by 
the carriers’ scattering with optical phonons, whereas 
e-e scattering is mainly responsible for the initial pho¬ 
tocarrier redistribution along the Dirac cone. The lat¬ 
ter observation is in stark contrast with the case of 
conventional direct-gap semiconductors, like GaAsl3l|. 
where the aniso trop ies are typically lost via carrier- 
carrier scattering [32|, l33|. 

Experiment — We study a chemical vapor deposited 
single layer graphene (slg)0,Is3 transferred onto a thin 
fused silica substrate, as described in Appendix. Struc¬ 
tural quality, uniformity and doping of SLG before and 
after transfer are investigated by Raman spectroscopy [35, 
Hi}. This shows that no damage occurs to the sam¬ 
ple as a consequence of the transfer process. The Ra¬ 
man measurements indicate that the sample is p doped, 
with a Fermi level~200meV[37l j38|. The D to G ra¬ 
tio measured at 514.5nm is ~0.24, which, combined 
with the estimated doping, corresponds to a small de¬ 
fect concentration^l0 11 cm~ 2 Idil. 

Transient absorption measurements are done with a 
15fs pump pulse at central photon energy of 1.62eV 
(A=765nm), see Appendix. A degenerate configuration is 
employed, where the same pulse is split and used for both 
excitation and probing. The experiments are performed 
as shown in FigQJ). The polarization of the probe pulse 
is rotated by 45 with respect to the pump. We sep¬ 
arate parallel and orthogonal components of the probe 
electric held with a polarizing beamsplitter (BS) after 
interaction with the sample. Our method ensures tem¬ 
poral synchronization of the different signals, as well as 
maximum spatial overlap for both measurements. This 
is crucial, since we target the quantitative comparison of 
the signals for both polarizations. 

The differential transmissions for parallel and perpen¬ 
dicular probe components are in Fig[2j Excited carriers 
result in photobleaching of the direct transitions, with 
an onset time as fast as the duration of the pump pulse. 
When probing the orthogonal polarization we only ob¬ 
serve a signal that, at early times, is significantly weaker, 
as compared to the configuration with parallel polariza¬ 
tions. We assign this difference to the anisotropy of the 
electron distribution in the Dirac cone, resulting from 
the pseudospin selective excitation probability. Already 
at short time delays<50fs the two signals start to con¬ 
verge. The two probe directions become indistinguish¬ 
able after~60fs, indicating that the carriers’ distribution 
becomes isotropic in the k x — k y plane, while the ther- 
malization along the cone continues. 

Importantly, the time required for the distributions 
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FIG. 2. Relative transmission change AT IT after optical ex¬ 
citation, probed with parallel (||) and orthogonal (_L) polar¬ 
ization. Representative measurements for (a) high and (b) 
low fluences. Insets: signals normalized to unity. 


to become isotropic does not depend on the excitation 
fluence. Pure e-e scattering would show a stronger de¬ 
pendence on the starting electron density, with increased 
probability of e-e collisions [i^. Indeed, e-e scattering 
across the Dirac cone occurs with a reduced phase space 
because each electron requires a companion with oppo¬ 
site momentum, as shown in FigJ3J:. This process is en¬ 
hanced at higher excitation densities when more scatter- 
ers are available. In contrast, the probability for phonon 
emission does not depend on carrier concentration, when 
the excitation density is far from saturation. The opti¬ 
cal phonons can scatter the electrons into any direction 
reducing their energy by hcoo and providing conserva¬ 
tion for any momentum, see Fig[3j:l, because the opti¬ 
cal branches show limited dispersion in the fc-space 411. 
This, combined with the characteristic timescale of a few 
tens of fs, allows us to identify e-ph scattering as the 
main isotropization process. 

To fully track the dynamics we perform a set of fluence- 
dependent measurements in which the excitation den¬ 
sity is varied between 14 and 705ytJ/cm 2 , see Fig|T] We 
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FIG. 3. Fundamental assumptions of our model, (a) Pump 
and probe pulses are approximated by rectangular profiles. 
They share the same spectrum and duration but differ¬ 
ent polarizations, described by the angles 9 1 and 62 with 
A 9 = 62 — 9i. The time delay between the two pulses is 
At. The model is valid in the case of non-overlapping pulses, 
(b) The photoexcited carriers relax towards the neutrality 
point by emitting optical phonons with energy Tllo 0 . The elec¬ 
tron cooling becomes inefficient at the bottleneck energy near 
huJo/2. E-e scattering does not change the total energy of the 
electronic subsystem, (c) Two-electron scattering across the 
Dirac cone shown schematically in cross-section at energy E. 
(d) Optical phonons scatter the electrons into any direction. 
A lower-lying cross section of the Dirac cone at energy E — Tno 0 
is shown by a thinner circle. 


evaluate the ratio between signals for parallel and per¬ 
pendicular polarization at different pump-probe delays 
At < lOOfs. On this time scale, we observe photobleach- 
ing, due to strongly nonequilibrium and hot (thermal- 
ized) carriers, with just the former being anisotropic. 
The ratio has a strong dependence on the carrier density 
for fluences<100/rJ/cm 2 , i.e. below 10 13 cm~ 2 carriers. 
Increasing excitation fluence, the isotropic contribution 
becomes dominant, and we observe a smaller difference 
between the signals arising for opposite polarizations. 

The experiments also show that the maximum differ¬ 
ential signal for perpendicular probe photons does not 
coincide with that for parallel. This depends on fluence 
and is assigned to the combination of scattering processes 
along the cone driven by e-e collisions, and momentum 
isotropization driven by phonons (25|. The photoexcited 
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FIG. 4. Ratio of relative transmission changes for parallel and 
perpendicular probe polarizations for a pump-probe delay of 
15fs, when the pulses are negligibly overlapping in time (At > 
St), but the occupation anisotropy is still substantial (At < 
r), see Fig[jI]for notations. Experiments (blue triangles) are 
in good agreement with Eq.Jl]) (red line). The inset is the 
experimental ratio at time delays between 5 and lOOfs. 
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carriers are first redistributed across the Dirac cone due 
to the optical phonon emission, Figj3jl, then along the 
cone in the new momentum direction, mainly via e-e 
scattering, Figj3}o. The experimental fingerprints of these 
phenomena are seen in Fig^insets), where the AT/To 
maxima for different probe polarization do not coincide, 
having a^lOfs delay. This helps once more to unveil the 
different role of e-e and e-ph scattering in the overall re¬ 
laxation dynamics. 


Our experiments also show that, at pump 
fluences<50^J/cm 2 , SLG has increased absorption, 
i.e. the transient signal becomes negative after~ lOOfs, 
as for Figl^b). Refs.[42l. |431] assigned this to phonon- 
assisted intraband absorption of the probe photons. Its 
contribution is stronger for lower excitation densities, 
where the weight of intraband transitions of a few hot 
electrons can be significant [ 42 ], |43| . At higher pump in¬ 
tensities, this is hidden under the dominating interband 
absorption involving a dense hot electron distribution. 


Model —Charge carriers in graphene near the K point 
can be described by the massless Dirac Hamiltonianj2*j| 
Hq = hvocr ■ k, where Kk is the two-component momen¬ 
tum operator, er is the pseudospin operator derived from 
the Pauli matrices, and vq = 1.05 • 10 6 ms -1 [I!]. The 
pseudospin orientation in the eigenstates of Ho depends 
on the direction of k. It is shown in Fig[T]as green arrows. 
The photon-carrier interaction is described by the Hamil¬ 
tonian Hi nt = ■ A, where A = A 0 cos (ujt — qz) is the 

vector potential created by the linearly polarized elec¬ 
tromagnetic wave E = Eo sin(u;t — qz) with ui = ‘I-kc/X 
the radiation frequency, and Eo = We assume nor¬ 
mal incidence q 1 k without momentum transfer from 
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photons to electrons. Note that the interaction Hamil¬ 
tonian H lnt inherits the pseudospin dependence from Hq 
and constitutes the interband transition selection rule of 

Fig HI 

We derive the photocarrier generation rate from the 
Liouville-von-Neumann equation [4^] for the density ma¬ 
trix written in the eigenstate basis of H 0 . The kinetic 
equation is solved within the duration of the pump pulse 
to < t < 0, and subsequently in the absence of exci¬ 
tation photons for t > 0 employing the relaxation time 
approximation and the solution at t = 0 as an initial 
condition, see Figj3li and Appendix. In order to quan¬ 
tify the contribution of the isotropic part of the pho¬ 
tocarrier distribution, we assume optical phonon emis¬ 
sion as the main mechanism responsible for decay of 
anisotropy within the first tens fs[25[. The electrons are 
relaxing to the bottleneck energy, see Fig j3}o, at which 
further cooling is strongly suppressed. Theory [-hi] and 
THz measurements [47| suggest that the e-h recombina¬ 
tion is a much slower process (ps time scale ji^, Irij'l. 
therefore, the photocarrier concentration is assumed con¬ 
stant. These approximations allow us to write the energy 
balance equation with only one free parameter, the hot 
electron temperature, and solve it analytically. 

The absorption coefficient A is defined as the ratio of 
absorbed to incident intensity. In the excited state we 
sum A t and Ah, which describe the absorption due to 
strongly nonequilibrium carriers and hot electrons, re¬ 
spectively, and Aq denotes the equilibrium value. These 
three quantities are given in Appendix. With negligible 
reflection j48l|. the optical transmissions read T = 1 — A 
and To = 1 — Ao, resulting in a differential transmission 
T ~ o T ° = — A ^ o 4 ° 1 ^ Q ■ The differential absorption is: 


A - An 


-(1 + 2 
+ tanh 


cos 2 A 9) 
hui 


2Aqv‘^ < ^>i At 

e t 


4 T h 
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( 1 ) 


2 

where Aq = tt a with a—j- the fine structure constant, 
u> the radiation frequency, 5u> the spectral width of the 
pulse, the pump fluence, At the probe time delay, 
St the pulse duration, and r ~ 30 fs the “orientational” 
relaxation time deduced from Figj2j where the two curves 
meet at a time delay ~50fs. The rate of e-e scattering 
would increase with <f>i and diminish r. The hot carrier 
temperature, Th, can be estimated as (see Appendix): 


Eqs. Hl2l are not valid at At = 0. EqQ]does not formally 
vanish for At r, because we do not take into account 
the energy dissipation at such a long time scale. Within 
the temporal constraints of the model, the description of 
processes occurring at At < St is limited. In particular, 
the slight delay in the maximum of the pump-probe sig¬ 
nal, as shown in the inset of Fig[3j cannot be reproduced 
by Eq. ©• This does not hinder the description of the 
isotropization dynamics, experimentally seen over 50fs. 

Eq. (jT]) predicts a relative differential transmission^ 
10 -4 -lCT 3 at pump fluences between 10 and 100^J/cm 2 , 
in agreement with the measurements in Figj2] We also 
evaluate the ratio between differential transmissions for 
parallel and perpendicular polarizations, and plot this 
as a function of pump fluence in Fig[4j The ratio de¬ 
creases at higher fluences because the relative contri¬ 
bution of the polarization-dependent term in Eq. (JT]) is 
strongly suppressed at higher Th, making the overall ex¬ 
pression less sensitive to AO. Physically, the higher exci¬ 
tation density delivers more heat to the isotropic Fermi 
sea, whose contribution suppressed the strongly nonequi¬ 
librium anisotropic component. 

Conclusions — We employed fluence-dependent and 
polarization-resolved optical pump-probe spectroscopy 
to resolve and explore the dominating relaxation mech¬ 
anism for photocarriers in graphene at ultra-short 
times scales. We found that optical phonon emission, 
rather than e-e scattering, is responsible for momen¬ 
tum isotropization and pseudospin relaxation in the non¬ 
equilibrium photocarrier occupation, while the initial 
photocarrier redistribution along the Dirac cone in a 
timescale of tens fs is due to e-e scattering. We pro¬ 
vided an analytical framework for the qualitative under¬ 
standing of the carrier dynamics. Ref. [28| suggested that 
the light generated anisotropic distribution of carriers in 
momentum space can be observed in electrical measure¬ 
ments despite their relaxation on ultra-fast time scales, 
as recently reported in Ref. [50 ] . Our model explains why 
it is possible: the continuous wave laser in Ref. [2^] re¬ 
sults in low photocarrier densities and strong anisotropy, 
allowing the pseudospin-polarized photocarriers to be de¬ 
tected in graphene pn-junctions[28|. The development of 
an analytical framework for the description of anisotropy 
dynamics in graphene paves the way for the qualitative 
design of novel photodetectors. 


T h 


tt 2 ah 2 VqUi 
6C(3 )u> 


( 2 ) 


with £ the Riemann zeta-function [ffl. The limit of 
—» oo is excluded in our model as we assume a system 
with linear response in $i. Th cannot be defined until 
a significant fraction of the strongly nonequilibrium car¬ 
riers relaxes to the hot Fermi-Dirac distribution[20], i.e. 
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FIG. 5. Noncollinear optical parametric amplifier output 
spectrum used in our pump-probe setup. 


APPENDIX 


for the density matrix p. It is convenient to 
work in the interaction picture [H3|. where the 
density matrix and light-carrier interaction are 
given by p 1 (t) = exp(j-H 0 t)p s (t) exp(-j-H 0 t), 
H( nt 0 ) = exp{j-H 0 t)Hf nt (t)exp(-j-H 0 t) with H 0 

the unperturbed Hamiltonian. The superscripts S and 
I stand for the Schrodinger and interaction picture, 
respectively. The Liouville-von Neumann equation can 
then be written as [45| : 

% = -\[HUt),p I {t)}. (3) 

Its solution can be obtained using an iteration procedure. 
First, Eq.® is integrated to give: 

t 

/(*) -/(*o) = ~ J ^Wnt(0>/(>')]• (4) 

to 


Sample preparation 

A 35pm. Cu foil is first annealed at~1000°C under a 
20sccm flow of H 2 for 30minutes, followed by 5sccm of 
CH 4 . The H 2 and CH 4 flows are kept constant for 30 
minutes, after which the chamber is left to cool down 
for~3 hours. SLG is then transferred onto 170/rm thick 
glass substrates by a wet etching method [H 0 . A sacri¬ 
ficial layer of polymethyl-methacrylate (PMMA) is spin 
coated on one side of the Cu foil. The Cu/SLG/PMMA 
stack is then left to float on the surface of a solution of 
ammonium persulfate (APS) in water. APS slowly etches 
Cu, leaving the PMMA+SLG membrane floating. The 
membrane is picked up using the target glass substrate 
and, after drying, the PMMA is removed with acetone. 


Ultrafast pump-probe setup 


Transient absorption spectroscopy is measured with a 
Yb:KGW regenerative amplifier system operating at a 
50kHz repetition rate. The laser drives a home-built non- 
collinear optical parametric amplifier (NOPA) which de¬ 
livers an output spectrum spanning over 0.35eV at a cen¬ 
tral photon energy of 1.62eV (A=765nm, see Fig® 1511. 
We compress the pulses to a temporal duration of 15fs by 
means of chirped dielectric mirrors [Hlj|. The probe pulses 
are detected with two photodiodes (for parallel and per¬ 
pendicular polarizations) followed by a lock-in amplifier. 


Photocarrier generation rate 


We 

from 


derive 

the 


the photocarrier generation rate 
Liouville-von Neumann equation fjjj 


Here, to has the meaning of initial time, when the in¬ 
teraction is on. Second, Eq.® is inserted back to the 
right-hand side of Eq. ® to give: 

f) f)I 7 

- E- = --[HUt),P I (t 0 )] 

t 

-£T J (?)]]■ (5) 

to 

It is possible to look for higher-order terms in H ( nt . Here, 
we restrict ourselves to the expansion up to the second 
order in H? nt , corresponding to the linear response in 
terms of the radiation fluence. Now we transform Eq. ((5]) 
back to the Schrodinger picture. The left-hand side of 
Eq.® then reads: 


dp 1 

~dt 




dp s 

~dt 


+ -[H 0 ,p s (t)]. 


( 6 ) 


The right-hand side of Eq.® consists of two terms: 


0 )] o)e* H ° At ] (7) 


and 

t 

j (8) 

to 

0 

[ dt"[H? nt (t), e& Hot " [Hg t (t +1"), p s (t + t")]e"^ ot "], 

to—t 

where t" = t' — t. A linearly polarized electromagnetic 
wave propagating in the ^-direction with wave vector q 
and frequency w can be described by the vector potential 
A = Aq cos (ujt — qz) which relates to the corresponding 
electric field as E = Eq sin (uit — qz), where Eq = ujAq/c. 
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In the case of normal incidence, q does not have any in¬ 
fluence on the in-plane carrier momenta k. The interac¬ 
tion Hamiltonian Hf nt {t) can be deduced from the tight- 
binding effective Hamiltonian [ 2 ^] and written as [30]: 


Htat(t) = hi n t (e^+e--*), 

where 

u _ ev 0 E 0 ( 0 e~ ie \ 
int “ 2w £ e ie 0 )' 


(9) 

( 10 ) 


with tan0 = Aq v /Aq x the polarization angle. Substitut¬ 
ing Eq.(j9D into Eos. (171811 . and neglecting the fast oscil¬ 
lating terms of the form e ±zwt , we obtain the following 
equation for the density matrix: 

^ + ~[H 0 ,p s {t)\=g{t ), (11) 

where 

0 

9{t) = ~2 / dt" {(e iut " +e~ iut ") 

to—t 

x[h int ,e^ Hot "[h int , p s {t + t")]e~^ Hot "}^ , (12) 

is the photogeneration rate. Eas. (11 11121) are the operator 
equations. To get the corresponding expression for the 
distribution function we rewrite Eos. (111I12D in the helic- 
ity basis, i.e. the eigenfunction basis of the unperturbed 
Hamiltonian Ho[13- The resulting distribution function 
is a 2 x 2 matrix. We retain only its dia gon al elements, 
relevant for interband optical transitions [45| . Thus, in¬ 
stead of Eq. ED we get: 

= S± [/ ± M1, ( 13 ) 

where g± [f±(t )] is given by: 

0 

g± [f±(t)} = -4 

h 

to-t 

x [f±(t +1") - f T {t + t")\ h 12 h 2 i} ■ (14) 

Here, f±{t) is the electron distribution function in either 
valence (subscript or conduction (subscript “+”) 

band, Tifl = 2hvok is the interband transition energy, 
and hij are the matrix elements of /ii n t given by: 

j _ evoEo f cos (<f> — 0) -isin(<^-0)^ 

ij ~ 2 w V*sin(<A- 0 ) cos(<f>- 0 ) )’ [ b) 

where tan <fr = k y /k x . The product hi 2 h 2 i can be thus 
written as: 

hi 2 h 21 = sin 2 {(j)-0). (16) 

Eq. (1161) reflects the polarization dependency in the pho¬ 
tocarrier generation rate given by Eq. ED. 


I dt "{ 


(u—Q) _i_ — it" (w— Q) 


Description of the pump pulse 


The pump pulse excites the carriers over the Fermi 
sea described by the Fermi-Dirac distribution f± \ The 
generation rate then derives from Eq. ED: 

(i) _ 2 h^h^i sin [(£ - t 0 )(wi - H)] 

9± ~ n 2 wi-n 

x (/£ 0) - /£ 0) ), t > to (17) 

where the pulse is switched on at to < 0, see FigEI 
The evolution of the non-equilibrium distribution func¬ 
tion within the pump pulse can be written as: 

Off 1 ) fW 

, t 0 <t < 0 (18) 

at t 


where r is the relaxation time, and reads: 


jCL) = - fi) 


2 h$h$ 


(/i 0) - /f) 


(19) 


1 + t 2 (oj 1 - fi) 2 h 2 (u)i - H) 2 
x (r(wi — ft) cos[(wi — fi)(t — to)] — sin[(u;i — fi)(t — to)]} 
to < t < 0. 


One can prove by direct substitution that f± \ given by 
Ea. (fl9P . satisfies Eq. ED- 

Once the pump pulse is switched off at t = 0, Eq. ED 
becomes: 


df¥_ 

dt 
F (i) 




Its solution is f±\t) = f±\ 0)e */ T , where f±\ 0) i 


t > 0 


( 20 ) 




is 


given by (1191) at t = 0. The non-equilibrium distribution 
function at t > 0 thus reads: 


f ( X ) = t(wi - H) 


(l), ( 1 ) 

12 21 (/£ 0) -4 0) ) ( 21 ) 


1 + T 2 (wi - H) 2 h 2 (uJi - H) 2 
x {r(w 1 — f l) cos[(wi — H)(5t] — sin[(wi — H)dt]} e ~ t ^ T , 
t > 0, 


where St plays the role of the pulse duration. The en¬ 
ergy relaxation time is assumed to be much longer than 
r. Therefore, the total distribution function at t > 0 rep- 
resents the sum of the hot Fermi-Dirac distribution fl j_ 
and the non-equilibrium addition, Eq. ED- 

The prefactor can be rewritten in terms of 

pump fluence = (cE 2 St)/(8 tt) (with E x the electro¬ 
magnetic wave amplitude), so that Eq. (fl6l) becomes: 


h W h W - 
n 12 n 21 ~ 


2TtaTiVQ<&i 

ujfSt 


sin 2 ((f> - 6»i), 


( 22 ) 


where a = e 2 /{Tic) is the fine structure constant, and 
tan = E\ y /E\ x . Since St is much longer than the typ¬ 
ical time scale oj/ 1 determined by the optical frequency, 
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we consider the limit of St —> oo. r must be set to infin¬ 
ity at the time scale w^ -1 , because t > St by definition. 
Using the formula p32j: 


lim 

St—too 


4sin 2 [(wi — tt)5t/ 2] 
(wi — fl) 2 St 


27T(5(wi — U) 


(23) 


we get the approximated expression for Ea. (l2Tl) : 


(i) 47r 2 aug$i , 2 

X (4°’ - /f) «(w i - n)e""\ 


t > 0. 


(24) 


Description of the probe pulse 


We now consider the optical absorption A = Ah + A T 
of the probe pulse governed by the hot electrons (index 
H) and strongly non-equilibrium short-living photocar¬ 
riers (index t) created by the pump pulse. We first cal¬ 
culate Aq in the absence of the pump pulse. The total 
number of the optical interband transitions Go within 
the pulse duration St can be evaluated by integrating the 
generation rate Ea. dTll) over the time t: 

Go = J dt J dt” cos [t"(ui 2 — fi)] 

ti ti~t 

X (/<?> - /<?>) hfMS ( 25 ) 

where is the equilibrium carrier distribution function 
at lattice temperature, and: 


(2) ( 2 ) _ 2-Kahv%$2 


u 12 n 21 ~ 


uj%St 


sin (</> - 0 2 ), 


(26) 


with <J> 2 , W 2 , 02 the probe pulse fluence, frequency, and 
polarization angle. Taking the integrals in Ea. (l25l) : 


G 

x 


o — 


27TQ:1Iq$2 
TlUn 


sm 2 ((f> - 0 2 ) 


4 sin 2 \{u) 2 — £l)5t/ 2] 
(u> 2 — Q)St 


(/i 0) - /. 


(27) 


We again exploit the fact that St » w/ -1 and utilize 
Eq. (l23l) for the transformation of the second line of 
Eq. (l27l) into the ^-distribution. The absorbed fluence 
for a given valley/spin channel can be then calculated as: 


This result agrees with previous measurements j48j . 

We now consider the probe pulse absorption due to 
the hot carriers created by the pump pulse. The hot 
carriers are described by the hot Fermi-Dirac distribution 
f±*\ which should now substitute f±^ in Eas. (l25l27l) . 
Since the chemical potential is much smaller than the 
excitation energy we set former to zero. The occupation 
difference can then be written as: 

f (o) _ f (0) = tanh ) (30) 

and the absorbed fluence in the presence of hot carriers 
for a given valley/spin channel as: 

_ 7ra , (hu>2 \ _ . 

= _tanh(^ —J$ 2 , (31) 

where Th is the hot carrier temperature estimated be¬ 
low. The corresponding optical absorption in the one- 
color pump-probe setup wi = u) 2 = u> is then given by: 

A h = 7ro;tanh \ ■ (32) 


Finally, we calculate the optical absorption due to the 
strongly non-equilibrium time-dependent carrier distri¬ 
bution after the pump pulse. The interband transition 
rate can be evaluated from Eq. (ED: 


u 

9 r(t) = J 2 J dt " C0S [t"( w 2 - U)] 


t\~t 


fW {t + t")-ft>(t + t") 


h [2) h {2) 
rt 12 ri 21 • 


(33) 


Here, the out-of-equilibrium distribution is given by 
either Eq. (T2T[) or (|21)) . depending on the approximation 
used. In what follows we employ the latter because we 
have already utilized a somewhat similar approximation 
to derive A 0 and Ah- Note that because of 

the SLG e-h symmetry. The number of optical interband 
transitions within the probe pulse then reads: 


G T =- 


4r/| 1} (0) 


h {2) h (2) 
ri 12 ri 21 


Ti (<jJ 2 — fl) 1 + (W2 — H) 2 T 2 


((w 2 - H)re- At/T 


+e <1 ' /r (sin [8t{w 2 — H)] — {w 2 — H)r cos [St{w 2 — H 


(34) 


$n 


' 

47T 2 


dflfl _ 7 tq; 

-^-2-Gq(H, (t>) = — $2 (28) 


Here, we assume /i 0 "* = 1 and = 0 at O = W 2 - Ho 
should include the spin/valley degeneracy: 


Aq = 7TQ:. 


(29) 


and the fluence absorbed within this process can be ob¬ 
tained by integrating G T over the whole fc-space, and by 
subsequent averaging over the pulse duration. The latter 
makes the fast oscillating terms in the second line of (1511) 
vanish. We thus get: 


$ T = 


7r 2 a 2 Ug$i4>2 

2hujiuJ2Suj 2 St 


[2 + cos (2A0)] e -At / r , 


(35) 

















where A 6 = 62 — Q\, and Su = W 2 — uq. In the one- 
color pump-probe setup uq = ui 2 = w, and plays 
the role of energy uncertainty (estimated as St -1 ). The 
resulting probe pulse absorption due to the strongly out- 
of-equilibrium carriers created by the pump pulse is: 


A T 


27r 2 a 2 i>g<I>i 

huj 2 Su} 2 6 t 


(1 + 2 cos 2 Ad) e~ At/T . 


(36) 


HOT TEMPERATURE CALCULATION 

To estimate the hot electron temperature as a func¬ 
tion of the pump pulse we assume that (i) the e-h re¬ 
combination process is much slower than r and, there¬ 
fore, the photocarrier concentration can be considered 
as a constant on a time scale of a few tens fs; (ii) the 
nonequilibrium photocarrier occupation relaxes towards 
the hot Fermi-Dirac distribution mostly due to the opti¬ 
cal phonon emission with the frequency to 0 , resulting in 
a characteristic electron energy ~ Tluiq /2 after thermal- 
ization, see Fig [3] The photoelectron concentration n p h 
for a given spin/valley channel can be written as[46j: 

f d2k A 1 ) /o 7 x 

Uph ~ J 4 +^+ ’ ( 37 ) 

whereas for photoholes we have [i^| : 

PP1 = /S(W“). (38) 

Here, is given by Eg. (El) . The energy balance equa¬ 
tions for a given spin/valley channel become: 

n ph-y^ + J ^/f=/ kfi H) i (39) 


p p ^ + J ^^ofc(i-/i 0) ) = j ^hv 0 k(l-f iH) ). 

(4 °) 

Here f± f ( f± f) ) is the Fermi-Dirac distribution at the 
lattice (hot) temperature. The integrals in Eas. (1371141)1) 
can be solved analytically in the case of intrinsic SLG, i.e. 
at zero doping. However, our samples are not intrinsic 
with a doping of the order of lOOmeV. To proceed, we 
assume that the chemical potential /x > 0 is higher than 
the lattice (room) temperature Tl, but lower than the 
hot electron temperature Th■ For electrons, we have: 


d 2 k {H ) 
4^0 kf + 


3 C(3 )T h 

47T H 2 Vq 


(41) 


d 2 k ,,(0) 

4 ^ 2 hv ° k f+ 


6ttTl 2 Vq 


(42) 


and the electron energy balance reads: 


■KOthJo&i 

8 uq 


3 C(3)T3 


6tT? 1 2 Vn 


47T h 2 Vn 


( 43 ) 


Here ( is the Riemann d-function|49 
temperature is given by: 


The hot electron 


Th = 


n 2 ahV qLUq® 

6C(3)wi 


_v_v 

9C(3)y 


(44) 


The holes at /.i < 0 can be considered in a similar way. 
Assuming typical values $1 = 10 — 100 /iJ/cm 2 , /j, of the 
order of 0.1 eV we find that Th ranges from 800 to 1800 
K. The /x-dependent term in Eg. (El) contributes weakly 
to Th and can be neglected. The hot electron and hole 
temperatures are equal within this approximation. Phys¬ 
ically, intrinsic electrons, while being at lattice tempera¬ 
ture, do not contribute much to the energy balance, even 
though their concentration might be high. Thus, we ar¬ 
rive at Eq. ©■ 
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